Depletion of the nuclear Fermi sea 
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The short-range and tensor components of the bare nucleon-nucleon interaction induce a sizeable 
depletion of low momenta in the ground state of a nuclear many-body system. The self-consistent 
Green's function method within the ladder approximation provides an ab-initio description of corre- 
lated nuclear systems that accounts properly for these effects. The momentum distribution predicted 
by this approach is analyzed in detail, with emphasis on the depletion of the lowest momentum 
state. The temperature, density, and nucleon asymmetry (isospin) dependence of the depletion of 
the Fermi sea is clarified. A connection is established between the momentum distribution and 
the time-ordered components of the self-energy, which allows for an improved interpretation of 
the results. The dependence on the underlying nucleon-nucleon interaction provides quantitative 
estimates of the importance of short-range and tensor correlations in nuclear systems. 
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INTRODUCTION 



Recent experiments at Jefferson Laboratory have clar- 
ified several aspects of the role of short-range and ten- 
sor correlations in determining the properties of nucle- 
ons in the nuclear medium. In an (e, e'p) experiment on 
"'^^C, an unambiguous signature of the presence of high- 
momentum nucleons was identified [1]. In the domain 
of missing energy and momentum probed by the experi- 
ment, the amount of single-particle (sp) strength identi- 
fied corresponded reasonably to scaled theoretical predic- 
tions of the self-consistent Green's function (SCGF) cal- 
culation for ^^0 [2] and correlated basis functions (CBF) 
calculations for nuclear matter in the local density ap- 
proximation [3] . The underlying mechanism for the pres- 
ence of these high-momentum components is associated 
with the strong repulsion that nucleons experience when 
they are in close proximity, generating a suppression of 
the relative wave function in coordinate space. In turn, 
this repulsion is required by the experimental nucleon- 
nucleon (NN) phase shifts. An additional mechanism 
that provides high-momentum components is the action 
of the NN tensor force, mediated to a large extent by the 
exchange of the pion. 

In even more demanding exclusive two-nucleon knock- 
out experiments on ^^O to the ground state of ^"^C [4, 5], 
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direct evidence for the presence of short-range proton- 
proton (pp) correlations was identified, in reasonable 
agreement with theoretical calculations [6, 7]. Definitive 
evidence for the importance of the nuclear tensor force in 
generating high-momentum components was presented in 
Ref. [8] , where the ratio of knocked-out proton-neutron 
(pn) to pp pairs from ^^C was found to be around 20. 
Theoretical relative momentum distributions exhibit a 
similar enhancement of pn over pp correlations due to the 
tensor force in the domain of momenta probed in this Jef- 
ferson Lab experiment [9]. Possible implications for the 
physics of neutron stars were discussed in Ref. [10]. 

Particle number conservation requires the high- 
momentum components to be accompanied by a corre- 
sponding depletion of the nuclear Fermi sea. A charac- 
teristic feature of this depletion in nuclear matter is its 
essentially momentum-independent character, except in 
the immediate vicinity of the Fermi surface [11]. A sur- 
vey of calculations [12-15] of the momentum distribution 
of nuclear matter at normal density in Ref. [11] demon- 
strated that, for all realistic NN interactions then avail- 
able, different many-body techniques consistently pre- 
dicted a depletion of the nuclear Fermi sea of a little over 
15%. About one-third of this depletion results from ten- 
sor correlations. In addition, the influence of three-body 
forces on the depletion appears to be rather insignifi- 
cant [12]. An (e,e'p) experiment on 2°*Pb at NIKHEF 
in a large domain of missing energy and a momentum 
range corresponding to the mean-field Fermi sea con- 
firmed that a global depletion between 15 and 20% of 
proton orbits below the Fermi energy explains all the 
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measured coincidence cross sections [16, 17]. To put these 
results in perspective it is useful to note that the deple- 
tion of liquid ''He for very small momenta is considerably 
larger and reaches about 50% [18], while the depletion of 
the electron Fermi sea in closed-shell atoms is essentially 
zero [19]. 

While there is uncertainty about the precise nature 
of the short-range part of the NN interaction, there is 
now evidence from recent lattice QCD calculations that 
the features of a strong repulsive short-range core emerge 
from first principles, particularly when the pion mass is 
reduced towards more realistic values [20] . The presence 
of short-range correlations is thus corroborated by QCD 
simulations and strongly suggests that fully microscopic 
nuclear many-body calculations should continue to ad- 
dress their consequences in detail. Recently developed 
realistic NN interactions, like CDBonn [21, 22] or the 
chiral interactions [23] have in general softer cores than 
older potentials [24, 25]. Even modern local potentials, 
such as Argonne vl8 (AvlS) [26], have a soft core com- 
pared to more traditional interactions. By assessing the 
different results for the AvlS and CDBonn interactions 
in the present paper, we plan to develop a measure of 
the remaining uncertainty of the role of short-range and 
tensor correlations in nuclei. 

A particularly well-suited technique for this goal is the 
SCGF method [27]. Within this approach, a fully self- 
consistent treatment of ladder diagrams for the interac- 
tion between particles that propagate with respect to 
a correlated (by short-range and tensor effects) ground 
state are accounted for. Recently, calculations includ- 
ing full off-shell effects have become available at finite 
temperature [28--31], providing thermodynamically con- 
sistent results that preserve sum rules, like e.g. conser- 
vation of the number of particles [32, 33]. Earlier appli- 
cations of this approach at T = involved a discretized 
version of the spectral distribution and have challenged 
the conventional interpretation of the nuclear-matter sat- 
uration problem [34]. Finite- temperature effects are ex- 
pected to smooth out the momentum distribution near 
the Fermi momentum, not unlike the influence of pairing 
correlations [25] . Thermal effects hardly affect the small 
and very high-momentum content of the ground state, 
however. This expectation is confirmed in general by the 
present work, although it requires some qualification in 
the case of low density or high temperatures, as we shall 
discuss in this work. 

The study of the asymmetry dependence of short-range 
and tensor correlations is motivated to a large extent by 
the (future) study of rare isotopes with large neutron ex- 
cess. Experimental studies employing heavy-ion knock- 
out reactions already suggest that the removal proba- 
bility for the minority species is strongly reduced com- 
pared to shell-model calculations, whereas the majority 
species exhibits properties that are essentially mean-field 
like [35]. Similar tendencies, but smaller in magnitude, 
have been obtained for the spectroscopic factors of pro- 
tons from a combined dispersive optical model analysis 



of a sequence of Ca isotopes [36, 37]. Fitting a huge 
collection of elastic proton scattering data while includ- 
ing results obtained from (e, e'p) reactions for quantities 
below the Fermi energy, the emerging complex optical 
potentials for protons exhibit a striking increase in sur- 
face absorption with increasing nucleon asymmetry. This 
translates into a qualitatively similar tendency in spec- 
troscopic factors as obtained from the zero-momentum 
occupation of protons calculated in the bulk for isospin- 
polarized nuclear matter [24]. A major purpose of the 
present work is to clarify the role of various variables, 
such as the temperature, the density or the choice of the 
NN interaction, in the determination of the momentum 
distribution for isospin-polarized systems. In this fash- 
ion, we can assess the importance of the short-range and 
tensor contributions to the bulk properties of asymmetric 
nuclei. 

Recent applications of the SCGF method to isospin im- 
balanced matter have generated predictions for the deple- 
tion of the proton and the neutron Fermi seas at fixed to- 
tal density as a function of nucleon asymmetry [24] . Some 
intriguing results were reported in this work that require 
a deeper understanding. First, an increasing diflference 
between the zero-momentum occupation of the neutron 
and proton Fermi sea was reported at a temperature 
corresponding to 10 MeV. Naively, this difference is ex- 
pected to be associated with the decreasing (increasing) 
importance of the nuclear tensor force for neutrons (pro- 
tons) with increasing asymmetry, since neutrons alone 
do not experience the strong ^Si-^Di tensor coupling, 
while protons will. There are, however, other physical ef- 
fects that can influence this difference and it is therefore 
useful to investigate whether the calculated predictions 
depend on the choice of the realistic interaction or on 
the temperature. Another feature that requires a better 
understanding is the opposite density dependence of the 
depletion for symmetric and pure neutron matter (PNM) 
reported in [24]. For symmetric nuclear matter (SNM), 
a decrease in the depletion of the zero-momentum state 
with increasing density is obtained, while the opposite 
result pertains for PNM in the case of the CDBonn in- 
teraction [21]. The result for SNM differs from earlier 
calculations, like those reported in Ref. [15] for a sep- 
arable version [38] of the Paris interaction [39], which 
predicted very little density dependence of this quantity. 
We propose to clarify this issue in the present work by 
analyzing more carefully the ingredients that determine 
the density dependence of the depletion, paying partic- 
ular attention to the importance of thermal effects and 
to the choice of the interaction. To provide a better in- 
sight into these dependences, we have relied and further 
clarified an approximate relation between the momentum 
distribution, above and below the Fermi surface, and the 
energy derivatives of the time-ordered components of the 
self-energy. These components differ from the retarded 
quantities that have been used in finite-temperature cal- 
culations. 

In Sec. II, we present the relevant ingredients of the 
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formalism of the finitc-tcmpcraturc implementation of 
the SCGF method. Results for SNM and PNM are pre- 
sented in Sec. Ill, where the role of the temperature, the 
density, and the choice of realistic interaction (softer vs. 
harder core) are clarified. A deeper insight into the den- 
sity and isospin-polarization dependence of the results is 
obtained by employing the relation between the momen- 
tum distribution and energy derivatives of the different 
time-ordered components of the nucleon self-energy. This 
topic is discussed in Sec. IV. Results for isospin-polarized 
matter are analyzed in Sec. V and further interpreted 
with the help of the relations discussed in Sec. IV. Fi- 
nally, conclusions are drawn in Sec. VI. 



II. SELF-CONSISTENT GREEN'S FUNCTIONS 
METHOD AT FINITE TEMPERATURE 

A. Single-particle propagators 

The results that are presented in the following have 
been obtained with the SCGF method at finite tempera- 
ture. In quantum statistical mechanics, the expectation 
value of any operator X is given by a grand-canonical 
average. 



(1) 



where (5 = 1/T is the inverse temperature and fi is the 
chemical potential of the system [32, 40, 41]. The many- 
body eigenstates of the system, |n), diagonalize simulta- 
neously the Hamiltonian and the particle number opera- 
tors, 



H\n) = En\n) , 
N\n) = 7V„|n) , 



(2) 
(3) 



and their Boltzmann sum (the denominator of Eq. (1)) 
defines the partition function of the system, Z. These 
eigenstates contain all the micro- and macroscopic infor- 
mation of the system, but an ab initio calculation in- 
cluding all of them is intractable. The description of the 
system can be substantially simplified, without loosing 
physical information, by considering particular kinds of 
excitations on top of the thermal bath. At the sp level, 
for instance, all the information of a homogeneous quan- 
tum many-body system is encoded in the time-ordered 
one-body propagator. 



iGT{k;t- t') = (t \ak{t)al{t') 



(4) 



given by the thermal average of a product of Heisenberg- 
picture operators associated with the addition and re- 
moval of sp excitations with momentum k between times 
t and t'. The time-ordering is implemented by the Wick 
operator, T, which arranges the operators in a chronolog- 
ical order and incorporates a (—1) factor according to 



the even (odd) nature of the corresponding permutation 
[40, 41]. The ordering of the time variable determines the 
analytical structure of Qr in the Fourier space associated 
with the time difference, t — t' . 
The correlation functions, 



(5) 
(6) 



are equal to the time-ordered propagator in the corre- 
sponding t > t' and t < t' domains. Because of the 
absence of time-ordering, Q-^ and Q'^ are analytical and 
well-defined functions at all energies. The Lehmann rep- 
resentation for , 



^G>{k■,t-t') = {a,{f:)al{t')^ 
ig<{k;t-t') = - (al{t')ak{t)) , 



g>{k,u;) = 27r^- 



-0{En-nN„) 

Z 



\{m\al\n)\'S{u;-Em + E„), (7) 

provides a physical interpretation for this function, as 

the probability of adding a particle with momentum k 
on top of the thermal bath and ending up in any possible 
final state, m, as long as the energy difference between 
the states coincides with ui. A similar interpretation ex- 
ists for Q'^ in terms of the removal of a particle from the 
bath. For a system in thermal equilibrium, the correla- 
tion functions are connected by detailed balance. 



g<{k,u;) = e^^^-'''>g>ik,uj). 



(8) 



yielding the so-called Kubo-Martin-Schwinger (KMS) re- 
lation [32]. The sum of g^ and g^ defines the sp spec- 
tral function of the system, A{k, oj), which is individually 
linked to either of them by the relations, 

g<{k,L,) = f{u;)A{k,u;), (9) 
g>{k,u;) = [l-f{u;)]A{k,u;), (10) 

with /(w) = [l + e'^^'^"'')] the Fermi-Dirac distribu- 
tion. 

The retarded propagator 

ign{k; t-t') = Q{t - t') (ak{t)al{t')) , (11) 



is related to the causal propagation of perturbations in 
the system, as enforced by the presence of the Heaviside 
function in relative time. Its spectral decomposition in 
frequency space is given uniquely in terms of the spectral 
function 



gn{k,uj) 



f 

J — < 



du;' Aik.J) 
27r a;+ — w' ' 



(12) 



where the notation w± = to ± irj, with rj infinitesimally 
small, has been introduced. Decomposing Eq. (12) in 
real and imaginary parts, 

lmgn{k,u;) = -lA{k,u;), (13) 



RegTi{k,u) 



r — 

J-oo 



lsn.gn{k,uj') 



(14) 
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a direct connection between ImQ-ji and the sp spectral 
function is found. Moreover, a dispersion relation {V de- 
notes a principal part integration) links the real part of 
Gtz to the energy dependence of its corresponding imag- 
inary part. 

Similarly, the spectral decomposition of the time- 
ordered propagator can be written as the sum of two 
terms 



(15) 



each of them depending both on the spectral function and 
an additional phase space factor. The real and imaginary 
parts of the up, 



lmg^{k,Lo)=^-f{u;)A{k,u;), 



— 

-oo 



dw' Im(7|(fc,cj') 



(16) 
(17) 



and the down components, 
1 



lmg^{k,u) = ~-[l- f{Lo)]A{k,uj), (18) 



/oo J / 
— 
-oo 



duj' Im^|(A;, w') 



UJ — UJ' 



(19) 



are linked by independent dispersion relations. It is then 

easy to show that the retarded propagator is related to 
the time-ordered components via the expressions: 

lTng^{k,u) = -f{io)lTngn{k,w), (20) 
lmg^{k,uj) = [I- f{uj)]lxngn{k,uj). (21) 

In the zero-temperature limit, the phase space factors 
in the dispersion integrals guarantee that the integration 
domains for the up and down components are discon- 
nected. The integration in goes up to the Fermi sur- 
face, Lu = fi, while the dispersion integral for starts 
at this point. Since the poles of the up (down) com- 
ponent are in the upper (lower) half of the complex en- 
ergy plane, one can associate the first (second) term to a 
"hole" ("particle") propagator. In contrast, the disper- 
sion integral associated with g-ji receives simultaneous 
contributions from above and below the Fermi surface. 
As we shall see in the following, this generates rather 
different energy dependences for the real parts of these 
three quantities. Let us also note that the knowledge of 
the spectral function is enough to compute the real and 
imaginary parts of all the different sp propagators [42]. 

The definition of the time-ordered propagator, Eq. (4), 
can be generalized to purely imaginary times with the 
proper treatment of the imaginary time-ordering [32]. 
The following KMS relation arises 



gTik,T) 



(22) 



which demonstrates that is a quasi-periodic function 
along the imaginary time axis. This suggests the follow- 
ing discrete Fourier representation 



gT{k,T) = 



1 



-0^ 



S^e-''-^g{k,z,)., 



(23) 



where z,y = 



Zu = '^'^^^ip'^ + correspond to fermionic Mat- 
subara frequencies [40-42]. An analytical continuation 
of the coefficients in the Fourier transform, gik^zj)^ to 
continuous complex values of energy can be uniquely de- 
fined [43]. The spectral decomposition of this function 



g{k,z) 



f 

J — < 



Auj' Aik.J) 
27r z — Lo' 



(24) 



is again uniquely given by the spectral function. Note 
that the same function, g, contains information on both 
the time-ordered propagator (when computed &i z = z^) 
and the retarded propagator (at z = a;+). In general, 
the sums over Matsubara frequencies can be transformed 
into energy integrals for the retarded components using 
complex analysis techniques. 

In momentum-frequency space, the Dyson equation for 
the sp propagator, 

z - Y,{k,z) 



2m 



g{k,z) = i. 



(25) 



reduces to an algebraic equation in terms of the self- 
energy, Ti^kjz). As a consequence of this equation, the 
self-energy inherits the analytical properties of the prop- 
agator. In particular, it can be decomposed in equiva- 
lent retarded, time-ordered or "less/greater than" con- 
tributions [42]. Taking the z — *■ oj+ limit in this last 
expression leads to a Dyson equation relating the re- 
tarded Green's functions and the retarded self-energy. 
Since these retarded components form a closed set of 
equations by themselves, they have been extensively dis- 
cussed, particularly in the nuclear context at finite tem- 
perature [29, 44, 45]. If one needs to distinguish between 
particle and hole contributions, however, time-ordered 
components become essential [44] . The decomposition in 
terms of up and down self-energies also allows for a more 
natural connection with zero-temperature SCGF calcu- 
lations [46, 47]. 



B. Ladder approximation 

The retarded self-energy fulfills the dispersion relation 

'"^ duj'lmErAk.oj') 



(26) 



The first term in this expression corresponds to an 
energy-independent Hartree-Fock contribution 

^HF{k) = J2 {kk'\V\kk')^ n{k') , (27) 



where we have introduced the NN potential (properly 
antisymmetrized) and the momentum distribution, n(fc). 
The latter includes correlation effects and can be com- 
puted from the sp spectral function 



n{k) = v j "^Aik, 



w)/(w) , 



(28) 
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where = 4 (i^ = 2) accounts for the degeneracy of nu- 
clear (neutron) matter. The energy-dependent, disper- 
sive contribution to tlic self-energy describes many-body 
processes that go beyond the basic mean-field approxi- 
mation. In the ladder approximation, one includes an 
infinite series of collisions between a particle and a series 
of particles and holes in the medium, thus accounting for 
the two-body scattering problem in the medium [40, 41]. 
The imaginary part of the retarded self-energy in this 
approximation reads 

/OO 1 / 

A{k' ,oj') {kk'\lmT{oj + uj'^)\kk') , 

(29) 

and is given in terms of the in-medium interaction, the 
sp spectral function and phase space factors, including 
a Bose-Einstein distribution, b{n) = [e-'^^^-^") - 
The T-matrix in the medium fulfills a Lippmann- 
Schwinger-like equation 

(feiA;2|T(0+)|fe3fe4)^ = {kik2\V\k3ki)^ 

+ {kik2\V\k5ke)A0ii{h,ke;n+)x 



{k5ke\T{n+)\ksk4)^ 



(30) 



where the retarded Q'^j is obtained from the product of 
two spectral functions and a phase space factor describ- 
ing the intermediate propagation of particle-particle and 
hole-hole pairs 



J — ( 



doj dw' 
2^ "27 



A{k,u)Aik',uj')x 
1 - fjco) - fju') 



properly accounting for Pauli blocking effects at finite 
temperature [41]. Once the real and imaginary parts of 
the self-energy are computed, one can feed back this in- 
formation into the spectral function via the Dyson equa- 
tion 



A{k,w) = 



-2ImSTC(/c,w) 



['^ - £ - ^e^nik, to)] + [ImEnik, u;)f 

(32) 



which can in turn be used to compute a new Qjj. By 
iterating this procedure until convergence, the SCGF re- 
sult within the ladder approximation is obtained. The 
importance of self-consistency arises from the fact that 
sum rules are preserved [48, 49] and also that thermody- 
namical consistency is fulfilled [50, 51]. In addition, self- 
consistency represents a democratic treatment of all the 
particles considered in the problem: the one for which the 
self-energy is calculated, as well as the ones it interacts 
with. In all the previous expressions, a grand-canonical 



ensemble has been assumed, so that there are two ex- 
ternal, fixed variables, the temperature and the chemical 
potential. The latter is not particularly well suited for 
in-medium studies, so it is customary to supplement the 
SCGF equations with an equation that fixes the total 
density of the system 



p = vYn{k,ij) , 



(33) 



where we have highlighted the dependence of the mo- 
mentum distribution on the chemical potential via the 
Fermi-Dirac factor of Eq. (28). 

The numerical solution of the SCGF equations is a de- 
manding task, due to two major difficulties. The first 
issue is related to the possibility of a pairing solution be- 
low a certain critical temperature, Tc, since nucleons tend 
to form Cooper pairs in certain density regimes [52, 53]. 
Although the SCGF formalism is capable of dealing with 
the description of the superfluid phase, we will only dis- 
cuss results for the normal phase, since the effects of 
pairing are concentrated in a narrow energy region (re- 
lated to the gap) around the chemical potential [25] . Any 
conclusion drawn for energies and momenta away from 
this region should therefore also be valid in the pairing 
regime and, most interestingly, in the zero-temperature 
limit. The second numerical difficulty implementing the 
SCGF method is related to the requirement to solve the 
coupled nonlinear equations in a wide range of momenta 
and energies. Since no quasi-particle assumption is made, 
complete off-shell propagation effects are included. In 
this case, different regions in the energy and momen- 
tum domains have to be sampled simultaneously, partic- 
ularly the quasi-particle peak and the high-momentum 
and energy components of the spectral functions. The 
smoothing of structures associated with finite tempera- 
ture is helpful in this direction, but the problem is still 
formidable. It has taken some time before complete, sat- 
isfactory numerical results have been obtained for realis- 
tic NN interactions. We refer the reader to Refs. [28-31] 
for further details on how the SCGF is implemented in 
practice. 

Once a self-consistent solution has been obtained, one 
has access to the spectral functions of a nuclear system 
at a given density and temperature. In addition to the 
sp properties, the spectral function also determines to a 
large extent the macroscopic properties of the system. 
The total energy, for instance, can be computed from 
the Galitskii-Migdal-Koltun sum rule [54, 55]. An ap- 
proximation to the entropy beyond the quasi-particle ap- 
proximation can also be obtained via the Luttinger-Ward 
formalism [51, 56]. All the remaining thermodynamical 
properties of the system are therefore accessible and in- 
clude properly the effect of short-range and tensor corre- 
lations. 

So far, we have discussed uniform nuclear systems with 
a single species: the nucleon, in the case of SNM, and the 

neutron, in the case of PNM. For a fixed total density, 
one can switch from one system to the other by modifying 
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the relative concentration of neutrons and protons. The 
asymmetry parameter 



Pn + Pp 

is a measure of the isospin imbalance. The SCGF 
method within the ladder approximation can be gen- 
eralized to the case of asymmetric nuclear matter and 
partially isospin-polarized systems can thus be ana- 
lyzed [24, 57, 58]. Since the approach accounts for both 
the short-range and tensor correlations associated with 
the underlying NN interaction, it can be used to gener- 
ate quantitative predictions for the importance of differ- 
ent types of correlations in isospin asymmetric systems. 
In this sense, it is a unique theoretical tool. Other theo- 
retical formalisms either cannot be generalized to isospin 
asymmetric systems (as is the case of the variational ap- 
proach [59]) or lack adequate consistency constraints (as 
in the case of Brueckner-Hartree-Fock theory [60] ) . 

In the following, we shall employ different microscopic 
NN interactions to quantify the uncertainty related to 
their short-range and tensor properties. Most of the cal- 
culations have been performed with the CDBonn [21] 
and the Avl8 [26] potentials. These phase-shift equiv- 
alent forces are representative of two subsets of realistic 
NN potentials. The first is a boson-exchange potential, 
with a soft short-range core and a somewhat small ten- 
sor component, as indicated by a relatively low D-state 
probability for the deuteron. The latter is a local poten- 
tial, with a harder short-range core and a more significant 
tensor coupling. Although these differences do not affect 
the two-body scattering observables (which are identical 
for any phase-shift equivalent interaction), they do in- 
fluence the in-medium properties due to their different 
off-shell behaviors [61]. To assess the importance of the 
different components of the nuclear force, we shall also 
discuss results obtained from the systematically simpli- 
fied family of NN interactions obtained by the Argonne 
group [62]. Each of these forces (from vlS to v4') has a 
successively simpler operatorial structure, with their pa- 
rameters refitted to describe the two-body system at their 
corresponding level of simplicity. Finally, the Reid93 po- 
tential, fitted to the Nijmegen partial wave analysis, has 
also been used as a benchmark [63] . The momentum dis- 
tributions obtained with these modern interactions will 
also be compared to results for older realistic interac- 
tions [11], that have somewhat stronger repulsive cores. 
In all cases, our present results have been obtained with 
partial waves up to J = 8 in the Hartree-Fock and J = 4 
in the dispersive contributions to the self-energy. 

III. MOMENTUM DISTRIBUTION OF 
ONE-COMPONENT NUCLEAR SYSTEMS 

We start the discussion of the SCGF results by high- 
lighting the effect of correlations on the momentum dis- 
tribution in SNM. Interaction-induced correlations have 



p=0.16 fm■^ T=5 MeV p=0.16 im^, T=10 MeV 




Momentum, k/kp Momentum, k/kp 



FIG. 1: (Color online) Correlated momentum distribution ob- 
tained with the Avl8 interaction at p = 0.16 fm~^ and two 
temperatures: T = 5 MeV (left panel) and T = 10 MeV 
(right panel) . The dashed line corresponds to the momentum 
distribution of the Free Fermi Gas under the same conditions. 



a distinctive signature in the momentum distribution, re- 
moving strength at momenta below the Fermi surface 
and shifting it to high momenta. In addition to the cor- 
relations induced by the NN force, however, the SCGF 
momentum distribution is sensitive to the temperature 
of the system. As a matter of fact, thermal correlations 
produce a similar signature on the momentum distribu- 
tion: low momenta are depopulated, while high momenta 
are thermally occupied. The amount of low-momentum 
depletion and high-momentum population, however, is 
significantly different for both types of correlations. One 
can distinguish between these two types of effects by com- 
paring the fully correlated SCGF prediction with the Free 
Fermi Gas (FFG) results, which is only sensitive to ther- 
mal correlations [64]. To this end, in Fig. 1 we show 
n{k) of SNM for Avl8 and the FFG for p = 0.16 frn'^ at 
T = 5 MeV (left panel) and 10 MeV (right panel). At the 
lower temperature, the Fermi-Dirac distribution of the 
FFG deviates very little from a step function, Q{kp — k), 
which describes the n{k) of the FFG at zero temperature. 
At this temperature and density, one can conclude that 
the thermal effects are small in the deep interior of the 
Fermi sea. Thermal effects, however, modulate n{k) close 
to the Fermi surface, which is no longer discontinuous at 
kp- Temperature is also responsible for the fact that, 
even in this very degenerate regime, the contribution of 
the states below fc^^- to the density sum-rule is only 86%. 
The remaining 14% corresponds to thermally populated 
states above kp- For the correlated case, n{k) is rather 
flat below kp and it presents a sizeable depletion (that we 
will later on characterize by ?t-(0)). With the comparison 
with the FFG, one can conclude that dynamical NN cor- 
relations are mainly responsible for the depletion. Due 
to these effects, the contribution to the density sum-rule 
of the states below kp is reduced to 75%, i.e. 10% of the 
strength is shifted to higher momenta due to NN correla- 
tions. Consequently, momentum components of the wave 
function above kp are generated. The high-momentum 
tail in n{k) at this density and temperature provides 52% 
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FIG. 2: (Color online) Density dependence of the occupation 
of the lowest momentum state at T = 5 MeV for Avl8 (solid 
lines), CDBonn (dotted lines) and the FFG (dashed lines). 
The left (right) panel corresponds to nuclear (neutron) matter 
results. Note the different vertical scale of the two panels. 



of the total kinetic energy per particle, to be compared 
to the much smaller 25% in the FFG. 

At a larger temperature (right panel), thermal effects 
do not only modify the FFG n{k) around the Fermi sur- 
face, but they also produce a depletion deep inside the 
Fermi sea. The lowest momentum state is depopulated 
by a few percent, n(0) = 0.97. The shape of n(k) dif- 
fers significantly from the step function, documenting the 
loss of degeneracy at this temperature and density, even 
though it is still far from the classical Boltzmann momen- 
tum distribution. The fraction of particles that occupy 
states below /cf in the FFG decreases to 73%, while they 
contribute only the 53% of the kinetic energy per parti- 
cle. At this temperature, the momentum distribution of 
the correlated system deep inside the Fermi sea exhibits 
a depletion, 1 — n(0) = 15%, which is approximately the 
sum of the depletion associated with the dynamic NN 
correlations ('^ 13%) plus the one coming from the ther- 
mal distribution of the FFG (3%). The high- momentum 
tail is however substantially more important in the cor- 
related case and the fraction of particles in states above 
the Fermi momentum increases to 44%, while the contri- 
bution to the kinetic energy amounts to 56%. 

The depletion of the momentum distribution below 
the Fermi surface can be taken as a measure of the im- 
portance of both thermal and dynamical correlations in 
many-body systems. In the degenerate regime, where 
thermal effects are unimportant, the amount of strength 
removed at low momenta is closely related to the struc- 
ture of the underlying NN force and is particularly sen- 
sitive to the short-range core as well as to the tensor 
components. In the correlated case, the low-momentum 
region of the momentum distribution is rather flat and 
almost independent of momentum as long as the temper- 
ature is relatively low and the density is large enough. 
For high temperatures or low densities, n(k) is softened 
by thermal correlations. One can therefore focus on the 
k = state to discuss generic depletion effects in both 
regimes. Figure 2 displays the density dependence of the 



occupation of this state for SNM (left panel) and PNM 
(right panel). Results for the Avl8 (CDBonn) interac- 
tion are shown in solid (dotted) lines, while dashed lines 
represent the FFG depletion under the same conditions. 

The non-interacting depletion is extremely helpful in 
understanding the influence of thermal effects on the 
low-momentum components. At high densities, degen- 
eracy dominates over temperature, and n{0) = 1 as ex- 
pected from the zero temperature FFG momentum dis- 
tribution. This degenerate regime determines the region 
where interaction-induced correlations represent the ma- 
jor contribution to the depletion in the correlated case. 
Consequently, the results obtained in this regime can be 
taken as faithful representatives of zero-temperature cal- 
culations. At lower densities, however, the FFG drops 
rapidly as density decreases, indicating that thermal ef- 
fects dominate the low-momentum components of n{k). 
The decrease observed in n(0) for the correlated case 
should therefore be associated with temperature effects 
rather than with NN correlations. In this region, the 
correlated n{0) yields interaction-independent results, as 
expected from the virial expansion [65] and illustrated 
in the figure. Note also that the onset density for the 
thermally dominated region is larger for symmetric mat- 
ter {p < 0.16 fm~^) than for neutron matter (p < 0.08 
fm~'^). This corresponds to the notion that, for a given 
density, neutron matter is more "degenerate" , i. e. has a 
larger Fermi momentum for the same value of p. 

If we take the deviation of n(0) from 1 as a measure 
of correlations, we can say that neutron matter is "less 
correlated" , in general, than nuclear matter (note the 
different vertical scale of the two panels of Fig. 2). More- 
over, the behavior of n{0) as a function of density in 
the region where thermal effects are unimportant is very 
different for PNM and SNM. In neutron matter, as the 
density increases, the population of low momenta drops. 
This result agrees with the intuitive idea that, as par- 
ticles are closer together on average, the effect of the 
short-range core increases and low-momentum strength 
is shifted to high momenta. Also in accordance with this 
picture, the depletion is more important for an interac- 
tion with a harder core (AvlB) than for a softer force 
(CDBonn) [25]. In stark contrast, for SNM the opposite 
density dependence is observed: the depletion is constant 
or decreases as the density increases. For AvlS, the de- 
pletion saturates at a constant value of n(0) ~ 0.87 as 
density increases, while for CDBonn a clear increase of 
n(0) with density is observed. 

What causes such a different behavior in the density 
dependence of n(0)? The major difference between neu- 
tron and nuclear matter lies in the role of the strong 
■^5*1 Di tensor component, which is only active in 
the latter. It appears that the effect of this component 
is twofold. On the one hand, it increases, in general, 
the depletion as compared to the effect of short-range 
correlations by themselves, as seen by comparing SNM 
and PNM results. On the other hand, tensor compo- 
nents seemingly modify the density dependence of n(0). 
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FIG. 3: (Color online) Momentum distribution of symmetric 
(upper panels) and neutron (lower panels) matter obtained 
with the Avl8 interaction at T = 5 MeV for different densi- 
ties. The right panels focus on the high-momentum region. 



This nontrivial result requires a deeper understanding for 
which the Green's functions formahsm offers a unique 
perspective by connecting directly the momentum dis- 
tribution and the in-medium self-energy, as will be dis- 
cussed in Sec. IV. Suffice it to say for now that the effect 
of the tensor force is not only limited to the depletion low- 
momentum components but is also responsible in higher 
order for binding effects, as illustrated by the deuteron. 

The normalization of the momentum distribution to 
the total density, Eq. (33), forces the missing strength at 
low momentum to be shifted to high-momentum states. 
The shift in strength may depend on the physical proper- 
ties of the system (temperature, density, isospin content). 
Figure 3 shows the momentum distribution of symmet- 
ric (upper panels) and neutron (lower) matter for Avl8. 
In all cases the temperature is 5 MeV and four different 
densities are displayed. In the left panels the momentum 
distribution is displayed as a function of the ratio k/kp, 
where kp = [Qir'^ /vp)^/^ is the Fermi momentum of the 
system. This allows for a clear inspection of the density 
dependence of n(fc) and can be instructive in understand- 
ing the scaling of the distribution as a function of density. 
At a constant temperature, an increase in density leads 
to a more degenerate system. Consequently, the high 
density results have a more zero-temperature-like struc- 
ture, with an almost constant depletion below kp and 
a substantial jump in n{k) at fc = kp. In contrast, at 
low densities, the momentum distribution is modulated 
by temperature at all momenta. Comparing the four dif- 
ferent densities, one can observe that one major density 
effect is the redistribution of strength in the region close 
to the Fermi surface. Note that on the present scale 



the effects on the depletion shown in Fig. 2 are hard to 
see. In general, though, Avl8 leads to a saturation of 
the depletion in symmetric matter and to a decrease of 
the population of low momenta for high-density neutron 
matter. 

The high-momentum components are shown in the 
right panels of Fig. 3 on a logarithmic scale for momenta 
ranging between 500 and 1500 MeV. In both symmet- 
ric and neutron matter, the high density results lead to 
a larger population of high-momentum states. This is 
in contrast to the low-momentum states, which exhibit a 
very different density dependence for the two cases. All in 
all, this suggests that the mechanism that produces high- 
momentum components only scales with density and is 
rather momentum independent. There have been sugges- 
tions in the literature that high-momentum components 
scale with fc|, as a function of density. Such results have 
been obtained for a separable version of the Paris inter- 
action [15] and a central (no tensor effects) version of the 
Reid soft-core interaction [66]. These interactions have 
a stronger core, which may explain the different behav- 
ior. In addition, in the limit of extremely dilute matter, 
the high-momentum components scale with k% for the 
hard-sphere Fermi gas [67]. Our calculations, based on 
more recent, softer realistic interactions, scale more like 
kl [30]. 



IV. RELATION OF n[k) TO THE 
TIME-ORDERED SELF-ENERGY 



The Green's functions method is unique in that it can 
provide an interpretation of the previous results in terms 
of the sp spectral function. The interplay between ther- 
mal and density effects in n(fc) can be qualitatively un- 
derstood in terms of the evolution in temperature and 
density of both the spectral function and the phase- 
space distribution (see Eq. (28)) [31]. Unfortunately, it 
is not straightforward to go beyond this sort of reasoning 
and find a simple connection to the in-medium interac- 
tion. A different perspective and additional insight can 
be achieved by considering an approximate relation be- 
tween the momentum distribution and the components 
of the time-ordered self-energy. The analysis will be per- 
formed for T = to make contact with earlier work on 
this subject [66, 68-71]. Let us write the Dyson equation 
at r = in the form 



(35) 



n=0 



where g^ is the time-ordered noninteracting sp propa- 
gator. The momentum distribution can then be written 
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n{k) = 



2Td 



(36) 



with C representing a contour composed of the real axis 
and closed by a semi-circle in the upper half-plane. In 
order to proceed it is necessary to decompose the self- 
energy into its different time-ordered contributions 



SffF(fc)- 



= Y.HF{k) + Sx(/c,tj) 
°°dw' Iml]x(fc,w') 



TT 



^^ik,Lu)^ (37) 
^^dw' Iml]|(fc,w') 



TT 



UJ' 



The first term in the expansion of Eq. (36), corresponding 
to Q'^ , yields the step function momentum distribution, 
n{kf = e{kF - k), of the FFG. Using Eq. (37), the sec- 
ond term (n = 1) in Eq. (36) yields three contributions, 
each of them with two noninteracting propagators, Q^. 
For k < kp, the contribution with the Hartree-Fock term 
vanishes on account of the double pole inside the contour, 
which leads to the derivative of a constant when Cauchy's 
integral formula is applied. The contribution with can 
be shown to be zero as well by using the corresponding 
dispersion relation for S| and applying the residue theo- 
rem (the residues cancel). Only the term with S| finally 
contributes and yields: 



n|(fc) = d^ReY,i{k,uj)\^=^o(^k) , k < kp ■, 



(38) 



by applying Cauchy's integral formula. The derivative 
is applied at the pole of the noninteracting propagator. 
Since the peak of the dressed propagator will normally 
occur at the quasi-particle energy 



e{k) 



2m 



ReS(fc,£(fc)), 



it is convenient to include an auxiliary potential 
U{k) = Rel](fc,£(fc)), 



(39) 



(40) 



anticipating that the derivative in Eq. (38) must be taken 
at the location where most of the sp strength is con- 
centrated. This leads to an additional constant term 
in Eq. (37), where U must be subtracted. An approx- 
imate result for n{k) for k < kp can now be obtained 
by expanding only up to first order at e(fc) in all 
contributions to higher-order terms {n > 1) in Eq. (36). 
This approximation will therefore be best for momenta 
as far as possible from the Fermi energy, since from the 
corresponding dispersion relation it is clear that a linear 
approximation to the energy dependence will be most ap- 
propriate for fc = (for which £(0) is farthest from e^). 
The series with only E^ terms yields 



n|(fc) = 



1 



1 - 9a;ReEx(A:,tj)|^^g(fe) ' 



k<kF, (41) 



Argonne vl8 



p=0.16fm , T=5 MeV 




0.5 1 1.5 
Momentum, k/k„ 



2 200 400 600 800 1000 1200 
Momentum, k [MeV] 



FIG. 4: (Color online) Correlated momentum distribution 
obtained with the Avl8 interaction at p = 0.16 fm"'' and 
T — 5 MeV (solid line). The ni{k) and n|(fc) contributions of 
Eqs.(41) and (43) are shown in dashed and dot-dashed lines, 
respectively. 



since a geometric series is produced with the energy 
derivative appearing in each higher order. Other terms 
including the constant T^hf and Ef will not contribute 
in a similar fashion as in the term with n = 1, as long as 
the validity of this expansion for E| is assumed. 

For momenta above kp, the noninteracting propaga- 
tors in Eq. (36) have poles outside the contour C and 
cannot contribute. Only terms with the derivative of S| 
will occur when a similar expansion as for fc < fc^? is em- 
ployed. The n = 1 term accordingly generates the first 
contribution to n(fc) for fc > kp^ given by 



nj(fc) = -ac^ReE|(fc, w)|j^=e(fc) , k > kf 



(42) 



For high momenta this will be the dominant contribution 
to n{k), since the location of e(fc) can be quite far from 
ep. Summing the geometric series from higher-order con- 
tributions with E| (with the approximation of keeping 
only linear terms in the energy) then yields 



niik) = 1- 



1 - ac^ReEt(fc,u;)|^=£(fc) ' 



k>kp. (43) 



While Eqs. (41) and (43) have not been rigorously derived 
here, we note that these results have been employed be- 
fore in the literature [70]. Indeed, it turns out that their 
numerical implementation generates an accurate approxi- 
mation to the full result for n(fc), except in the immediate 
vicinity of the Fermi surface, as will be illustrated in the 
following. Moreover, these equations allow for a deeper 
understanding of the observed trends in the depletion of 
the Fermi sea. Since the decomposition of the self-energy 
in terms of "up" and "down" components is essential for 
this derivation, it clarifies the need for introducing these 
terms also at finite temperature. 

For a self-consistent calculation at finite temperature, 
expressions (41) and (43) reproduce the momentum dis- 
tribution obtained from the complete spectral function 
extremely well in their respective momentum range. An 
example is shown in Fig. 4 for Avl8 at p = 0.16 fm~'^ and 
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FIG. 5: (Color online) Different components to the imaginary 
(upper panel) and real (lower panel) parts of the self-energy 
for CDBonn at p = 0.16 fm"^ and T = 5 MeV. 



T — 5 MeV. The hole contribution coming from Eq. (41) 
(dashed hne, left panel) is a very good approximation to 
the momentum distribution for momenta below 0.75kp. 
Similarly, for momenta above 1.25kp, n{k) is very well 
approximated by Eq. (43) (dash-dotted line, right and 
left panels). A similar agreement is found in a wide 
range of densities and temperatures, and also for other 
NN forces. We note that, for these results to be valid, it is 
essential to calculate the derivatives at the corresponding 
quasi-particle energies and therefore self-consistency is a 
necessary ingredient in obtaining this good agreement. 

Having found a direct connection between n[k) and 
the self-energy, one can try to establish a relation be- 
tween the density dependences of n{k) and the density 
and energy dependence of the self-energy. To this end, 
let us discuss the structure of and Figure 5 shows 
the self-energy of the k — state as a function of energy 
for the CDBonn potential at p = 0.16 fm~'^ and T = 5 
MeV. The upper panel displays three different compo- 
nents of S. The solid line corresponds to the retarded 
self-energy, which has traditionally been used and dis- 
cussed in finite-temperature SCGF calculations [30, 44]. 
Note that ImE^j, is negative at all energies and displays 
two distinctive peaks: a relatively small peak in the 
"hole" region, uj < fj,, and a much larger structure in 
the "particle" region, uo > ji. While the hole components 
decay rapidly for negative energies, the particle side has 
a much slower decay in energy due to the combination 
of the large particle-particle phase-space and short-range 
correlation effects, in complete analogy with earlier re- 
sults obtained for T = [46]. 

Since the self-energy has the same analytical structure 
as the sp propagator, the time-ordered contributions of 
E can be obtained from the retarded component by using 
equivalent expressions to Eqs. (20) and (21). The phase 



space factors in the imaginary parts isolate the respec- 
tive particle and hole regions. Consequently, ImS| and 
ImEj have a single peak below and above the Fermi sur- 
face, respectively. In the case of the hole contribution, 
the peak has the opposite sign compared to the retarded 
component, but for the particle one it coincides with the 
ImE^j, peak. The differences between the retarded and 
the corresponding time-ordered imaginary components 
will be larger in the lo ^ fj, region, due to the variation 
of the thermal distribution in this region. In particular, 
at higher temperatures, the opening of phase-space can 
eventually lead to the appearance of additional peaks in 
the self-energy of each component, due to the overlap be- 
tween the "opposite" peak and the non-zero Fermi-Dirac 
distribution. 

The different imaginary components of E are related to 
their respective real parts by the dispersion relations of 
Eqs. (14), (17) and (19). These real parts are illustrated 
in the lower panel of Fig. 5. The important consequences 
of the isolation of the hole and particle components in 
the imaginary part can be clearly observed now. ReE| 
decays rapidly at low and high energies, because the hole 
strength in the imaginary part is limited to a relatively 
narrow energy domain. This component exhibits a wiggle 
at small negative energies, associated with the hole peak. 
The down component of E has a substantial contribution 
at all energies, as a consequence of the extremely long 
energy tail of the imaginary part [25, 30, 31]. It also 
has a peak structure, related to the onset of a very fast 
decrease of ImEj near uj ^ ^. This threshold, together 
with the relatively large and constant energy tail, are at 
the origin of a bell-shaped minimum. As we shall see 
in the following, it is this minimum that dominates the 
low-energy components of the momentum distribution. 
The real part of the down component has a zero at high 
energies, ~ 1.5 GeV, and it is still substantial at very high 
energies. These results are in good qualitative agreement 
with the zero-temperature calculations of Ref. [46] for a 
much harder interaction. 

The decomposition of E in its time-ordered compo- 
nents can also be used to better understand the struc- 
ture of the retarded self-energy. The sum of ReE^ and 
ReE| must coincide with the real part of the retarded 
self-energy [42]. The high positive and negative energy 
tails of the retarded self-energy have to be related to the 
down component of E, which is the only active contribu- 
tion in these regions. Moreover, the double-peak struc- 
ture of ReE-Tj at w ^ is caused by the superposition 
of the up and down peaks and wiggles in this region. In 
particular, the uj < fj, peak in RcE-r^ is a reflection of 
the negative side of the wiggle in ReEj , while the lo > ji 
peak is mostly due to the down component. The rel- 
ative importance of the hole and particle peaks changes 
with momentum, with the "particle" side becoming more 
prominent as k increases. In spite of this important shift, 
the connection between the up and down components and 
the high energy and peak structures of the retarded self- 
energy remains valid. 
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FIG. 6: (Color online) Imaginary part of E| as a function of 
energy for the A; = state. Different densities are displayed 
with different line styles. The upper panels have been ob- 
tained with CDBonn, while the lower panels correspond to 
Avl8. Symmetric (neutron) matter results are shown in the 
left (right) panels. Note the difference in scales between the 
upper and lower panels. 



The self-energy is related, via Eq. (41), to the mo- 
mentum distribution below kp. This suggests that the 
density dependence of the low-momentum components of 
n{k) can be understood in terms of the self-energy. Let us 
focus for simplicity on the fc = state. Figure 6 displays 
the energy dependence of ImSj for different densities in 
the case of symmetric (left panels) and neutron matter 
(right panels). The upper (lower) panels correspond to 
the CDBonn (Avl8) interactions. The imaginary part of 
Til has a very large structure in the uj > fi region that 
only decays at extremely high energies. The area covered 
by this function depends substantially on the density, the 
NN potential and the isospin content of the system. In 
the high energy regime, ImE^ displays a linear depen- 
dence in density for both nuclear and neutron matter. 

A significant difference is found by comparing upper 
and lower panels in Fig. 6: the amount of imaginary 
part shifted to high energies is substantially larger for 
Avl8 than for CDBonn (note the very different verti- 
cal scales). The down component of the self-energy is 
therefore very sensitive to the structure of the underly- 
ing NN force. This result suggests that hard interac- 
tions produce a much larger imaginary part than soft 
ones do [25], for both SNM and PNM. For a given NN 
interaction, however, the peak of ImE^ at high energies 
is relatively similar for the two systems, as seen by com- 
paring left and right panels. While in PNM the peak 
can be attributed basically to short-range correlations, 
SNM has an additional contribution coming from the ten- 
sor components that generates an important imaginary 
part below 1000 MeV, as discussed in Refs. [11, 72]. It 
is good to remember that the presence of an imaginary 
part in the self-energy above the chemical potential is 



FIG. 7: (Color online) The same as Fig. 6 for the real part 
of E|. The dots represent the position of the quasi-particle 
peak. 



directly responsible for the appearance of sp strength at 
those energies for momenta below kp- In nuclear mat- 
ter, the additional '^S'l-^Di coupled channel with both 
its short-range and tensor contribution is thus respon- 
sible for the stronger and more pronounced imaginary 
part starting above the chemical potential as compared 
to neutron matter, where it is absent. The position of the 
peak in ImE| also changes when going from symmetric 
to neutron matter, but in an interaction-dependent way. 
While for CDBonn the peak is shifted to high energies 
for neutron matter, the contrary happens for Avl8 [25]. 

The energy dependence of the imaginary part has an 
immediate impact on the real part of Ej^ by means of the 
dispersion relation. Figure 7 shows this component in a 
narrower energy range for various densities. In general, 
the real part is proportional to the density, with higher 
densities leading to more attractive real self-energies (in 
the range displayed in the figure). Comparing left and 
right panels, one observes an important difference be- 
tween symmetric and neutron matter. While a distinc- 
tive inverted peak develops in nuclear matter for the two 
NN interactions, in neutron matter only a very broad 
minimum is observed. The origin of this peak is associ- 
ated with the rapid drop, from zero to negative values, 
of ImE| in a narrow energy band around w ~ /i. This 
sharp threshold occurs for symmetric matter (see the up- 
per left panel of Fig. 6, for instance), but it is not as sharp 
in neutron matter, reflecting the important on/off role of 
the tensor force in these systems. As a consequence of 
these large differences, the dispersive counterparts of the 
threshold of the imaginary part are very different for both 
systems. 

The presence of the inverted peak in symmetric mat- 
ter has important consequences for the momentum dis- 
tribution at low momenta. To compute n{k), the partial 
derivative with respect to the energy of ReE| is computed 
at the quasi-particle peak. The quasi-particle energy for 
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FIG. 8: (Color online) Partial derivative of E| as a function of 
energy for the fc = state at T = 5 MeV. Different densities 
are displayed with different line styles. The upper panels show 
symmetric matter results, while the lower panels correspond 
to neutron matter. Left (right) panels have been obtained 
with the Avl8 (CDBonn) interaction. 



= at different densities is identified by dots in Fig. 7. 
As expected, at zero momentum the quasi-particles are 
more bound for higher densities, so their energies shift 
to more negative values with respect to the chemical po- 
tential. In SNM, this shift moves the quasi-particle con- 
tribution further away from the inverted peak, which is 
displaced to higher energies as density increases. The 
combination of the two effects (attractive quasi-particle 
shift, repulsive peak shift) is such that the slope at the 
quasi-particle energy is reduced with increasing density. 
For PNM, the quasi-particle shift is larger, but the in- 
verted peak is absent and the broad structure is such 
that the local slope around the quasi-particle contribu- 
tion decreases with density. 

To clarify these issues, we show in Fig. 8 the partial 
derivative with respect to the energy of the self-energy 
at T = 5 MeV for different densities. Upper (lower) pan- 
els correspond to symmetric (neutron) matter, while left 
(right) panels show results for Avl8 (CDBonn). Let us 
first focus on the SNM case (upper panels). At all den- 
sities, the inverted peak in the real part (a zero in the 
derivative) lies inside the uj > fj, region, and it is there- 
fore not visible in the figure. Instead, we have chosen to 
display the oj < jj, regime, where the derivative should 
be taken according to Eq. (41). This corresponds to the 
wings of the peak on the left side. At these energies, 
the slope of ReEj decreases with energy and, in gen- 
eral, has more negative values at higher densities. The 
increase of the slope with density also appears to be re- 
lated to the hardness of the interaction. While for CD- 
Bonn the slope changes very little as density increases. 
Avis leads to much larger variations in SNM. The soft 
interaction, therefore, generates a small variation in den- 
sity of 9[jReS| which, together with the quasi-particle 
shift, produces a decreasing depletion as the density in- 



creases. This observation therefore explains the counter- 
intuitive behavior observed in the left panel of Fig. 2. For 
a harder interaction, such as AvlB, the effect of density 
in the slope is much larger, while the quasi-particle shift 
is similar to CDBonn. Consequently, there is the possi- 
bility for the depletion to saturate (or even decrease) at 
high densities. Earlier results for a separable version of 
the Paris potential exhibit little density dependence for 
the depletion [15] not unlike our results for AvlB in the 
range from once to twice normal density. Note that the 
density dependence of the self-energy is closely related to 
the strong correlation effects induced by the short-range 
components. 

The density dependence of the partial derivative is also 
important in neutron matter (lower panels). Again, the 
slope changes much more as a function of density for 
Avis than for CDBonn. Due to the quasi-particle shift, 
the derivative at zero momentum decreases. Thus, the 
depletion increases with density, in agreement with the 
right panel of Fig. 2. For a given NN interaction, the 
overall decrease of cJ^^ReSi with density and the quasi- 
particle shift are qualitatively similar for both neutron 
and nuclear matter. The origin of the differences between 
SNM and PNM should thus be attributed to the energy 
dependence of the slope, which is much more pronounced 
for symmetric matter than for neutron matter. This is 
particularly true for energies close to a; ~ /i, where a 
sharp decrease is observed associated with the inverted 
peak in the real self-energy. As we have already dis- 
cussed, this is largely related to the threshold behavior 
of the imaginary part of Ej. In turn, this behavior must 
be partly ascribed to the tensor component of the one- 
pion-exchange interaction. 



V. MOMENTUM DISTRIBUTION OF 
TWO-COMPONENT NUCLEAR SYSTEMS 

Up to now, we have considered the momentum distri- 
butions for two extreme cases, namely SNM and PNM. 
It is also illustrative to study the dependence of the pro- 
ton and neutron momentum distributions, nr{k), as a 
function of the asymmetry parameter, a. This analysis 
should shed some light on the the isospin dependence of 
single-particle nuclear properties coming from volume ef- 
fects. In particular, tensor effects can be highlighted in 
asymmetric systems and can be pinpointed by a compar- 
ison of different NN potentials. 

The asymmetry dependence of the neutron and proton 
momentum distribution is presented in Fig. 9 at p = 0.16 
fm^'^ and T — 5 MeV. The left (right) panel contains the 
results for FFG (AvlS). For a = (symmetric case), 
the momentum distributions of neutrons and protons co- 
incide in both cases. Let us discuss first the left panel, 
where we show the momentum distribution of the isospin 
asymmetric FFG. For the symmetric case {a = 0) and 
as discussed previously in the context of Fig. 1, the ther- 
mal effects on n{k) are only visible at the Fermi sur- 
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FIG. 9: (Color online) Isospin asymmetry dependence of the 
Free Fermi gas (left panel) and the correlated momentum dis- 
tribution from Avl8 (right panel). The density (p = 0.16 
fm~'^) and temperature (T = 5 MeV) are fixed and the re- 
sults for different asymmetries are shown in different line- 
styles. Symmetric matter results correspond to the central 
solid curve, while curves to its right (left) correspond to 
neutron (proton) distributions for the corresponding isospin 
asymmetric cases. 



face and almost no depletion occurs inside the Fermi sea. 
The same is true for the most abundant component as 
asymmetry increases: the depletion is negligible and the 
Fermi momentum is displaced to larger values. For pro- 
tons, however, n(k) is clearly affected by temperature, 
exhibiting a large change in shape in a wide region of 
momenta. At large asymmetries, in particular, the mo- 
mentum distribution of the less abundant components 
differs substantially from the step function, loosing its de- 
generacy character and decreasing the occupation at the 
origin. In some sense, when increasing the asymmetry, 
the less abundant component moves closer to a classical 
momentum distribution, while the most abundant com- 
ponent becomes more degenerate. Note however that, at 
a — 0.8, one is still far away from the classical regime 
even though thermal effects are important. This anal- 
ysis suggests that one should take into account thermal 
corrections when analyzing the change in asymmetry of 
nr{k) also for the correlated case. 

Focusing on the right panel, one observes that the most 
abundant component (neutrons for positive a) gets less 
depleted when the asymmetry increases, i.e. neutrons 
become "less correlated" . This is in contrast to the FFG 
results, which, for neutrons, exhibit no change inside the 
Fermi sea. The decrease of depletion for the most abun- 
dant component in this density and temperature should 
be taken as a pure NN correlation effect. This behavior 
can be explained in simple terms as follows. Although 
the total number of pairs is the same when increasing 
the asymmetry at constant density, some of the pn pairs 
are replaced by neutron-neutron pairs. The latter corre- 
lations are weaker than the pn ones, due to the absence 
of tensor effects, and therefore neutrons become less cor- 
related at large a's. Conversely, the momentum distri- 
butions of the less abundant species (protons) become 
more depleted with asymmetry. A single proton sees an 



FIG. 10: (Color online) Isospin asymmetry dependence of the 
neutron (filled circles) and proton (squares) lowest momentum 
occupation. Correlated results for the Avl8 (solid lines) and 
CDBonn (dotted lines) interactions are compared to the Free 
Fermi Gas (dashed lines) predictions. 



increasing number of neutrons when the asymmetry in- 
creases, i.e. pp pairs are replaced by the more correlated 
pn pairs, which results into a more depleted proton mo- 
mentum distributions. Together with this, protons "feel" 
more the effect of temperature. The proton density in the 
asymmetric system decreases as a increases. At finite 
temperature, this is translated into the fact that np{k) 
becomes dominated by thermal effects at large asymme- 
tries and looks less degenerate. The plateau associated 
with particle states, for instance, which is clearly ob- 
served for neutrons at all asymmetries, gets washed out 
for protons at large a. Consequently, one should be care- 
ful in associating the increase of the depletion of protons 
only to dynamical NN correlation, because the momen- 
tum distribution of the less abundant component might 
be strongly influenced by thermal effects. 

The information about the isospin dependence of the 
depletion is summarized in Fig. 10, where nr(0) is plot- 
ted as a function of the asymmetry at T = 5 MeV for 
p = 0.16 fm-3 (left panel) and p = 0.32 fm'^ (right 
panel). The FFG (dashed lines) results are compared to 
the Avl8 (sohd hues) and CDBonn (dotted lines) SCGF 
predictions. Again, the FFG results provide a measure of 
thermal effects. On the one hand, the occupation of the 
most abundant component at zero momentum does not 
change when the asymmetry increases, indicating that 
this species is totally degenerate in the whole range of 
asymmetries. On the other hand, the corresponding oc- 
cupation of protons (the less abundant component) is 
close to 1 at small asymmetries and decreases as the 
asymmetry increases. At p = 0.16 fm~'^, this accounts 
for a - 2% effect for a ^ 0.6 and a - 11% for a. = 0.8. As 
a matter of fact, in the limit a — > 1, the protons become 
an impurity gas in a Fermi sea of neutrons, thus behav- 
ing as a classical gas with np(0) 0. The steepness of 
the change in np(0) depends on the total density and the 
temperature of the system. At a higher density (right 
panel) , the thermal effects affecting protons do not set in 
unless the asymmetry is very large and, even at a = 0.8, 
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the effect is only of ~ 2%. Let us stress the fact that 
this is a pure thermal effect: for the FFG at zero tem- 
perature, n-riO) = 1 for both protons and neutrons at all 
asymmetries. 

For the correlated depletion, the occupation of the zero 
momentum state is an increasing (decreasing) function 
of the asymmetry for the most (less) abundant compo- 
nent. The behavior is very similar for both NN interac- 
tions, although the occupations for Avl8 are systemat- 
ically smaller than those for CDBonn at both densities. 
The differences between the two potentials increase with 
density (see right panel) . The comparison with the FFG 
allows us to identify the regime in which NN correla- 
tions dominate over thermal effects. At T = 5 MeV and 
p = 0.16 fm^^, for instance, and up to asymmetries of 
about a ~ 0.4 — 0.6, the occupations of neutrons and 
protons can be attributed to dynamical correlations and 
should provide a good estimate for the case of zero tem- 
perature. Up to these asymmetries, both the neutron and 
the proton depletion change almost linearly with asym- 
metry. In the right panel, the depletion of protons starts 
to bend down for larger asymmetries, attributable to the 
onset of thermal effects. In a more degenerate case (right 
panel), where thermal effects are almost negligible, the 
asymmetry dependence of both tIt- (0)'s is again found to 
be linear. 

Asymmetries of stable nuclei belong to the lower range 
of asymmetries, i.e. a ~ 0.2 for ^"^Pb. The SCGF pre- 
dictions should be valid in this range, where thermal ef- 
fects are unimportant. Unfortunately, the difference be- 
tween the occupation of protons at this asymmetry and 
the symmetric case is only ^2%, too small to allow ex- 
perimental verification. Nuclei at larger asymmetry pro- 
duced at future rare isotope facilities may provide a bet- 
ter testing ground. Moreover, integrated effects over the 
whole neutron and proton density profiles might also en- 
hance the effect of asymmetry. Let us note that, for finite 
nuclei, one should also take into account surface proper- 
ties and their isospin dependence. In particular, an asym- 
metry dependent effect associated with surface properties 
has been identified for protons in Ca isotopes [36, 37]. 

A surprising feature arises when comparing the two 
panels of Fig. 10. In spite of the differences observed in 
both SNM and PNM for the Avl8 and CDBonn results, 
the asymmetry dependence of nr{0) is very similar at all 
densities. This is shocking because, as we have stressed 
so far, both forces have a rather different short-range be- 
havior and tensor structure. Given the almost linear de- 
pendence of nr(0) with asymmetry, a better insight into 
these differences can be gained by plotting the difference 
n-niO) — n-p(O). In Fig. 11, we present this "iso-depletion" 
for different interactions as a function of the asymmetry 
at T = 5 MeV for the same densities as before. The 
left panel shows that this difference is the same for a 
wide variety of modern NN potentials, independently of 
their short-range or operatorial structure. This seems 
to suggests that the iso-depletion is fixed by the phase- 
shifts, most probably via their isospin dependence. This 
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FIG. 11: (Color online) Difference of neutron and proton 
occupation of the lowest momentum state as a function of 
isospin asymmetry for different interactions. The two pan- 
els show results at T = 5 MeV for p = 0.16 fm"^ (left) and 
p = 0.32 fm"^ (right). 
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TABLE I: Occupation of the lowest momentum state for dif- 
ferent NN interactions in symmetric and neutron matter. 



is corroborated at higher densities by comparing the very 
similar Avl8 and CDBonn predictions (left panel). 

In both panels, there are two results that fall below the 
main iso-depletion line. Immediately below most of the 
interactions, one finds the results corresponding to the 
Av4' potential. This has an extremely simplified operato- 
rial structure, with only a spin-isospin part and no tensor 
components, and is fitted to reproduce the binding energy 
of the deuteron [62] . The fact that it lies significantly be- 
low the other results shows the importance of tensor ef- 
fects for isospin asymmetric systems. It appears that the 
tensor force tend to increase the difference between neu- 
tron and proton momentum distributions as asymmetry 
increases. Moreover, the comparison with Av6' and AvB' 
suggests that, once the tensor components are included 
in a force, the iso-depletion will remain almost the same 
independently of the extra spin-orbit terms. Let us again 
stress the fact that such an agreement for different poten- 
tials is very surprising, particularly if one considers the 
fact that the momentum distributions of neutrons and 
protons can be different for each potential. 

The FFG iso-depletion is, in all the cases explored, be- 
low the correlated results, which shows the importance 
of NN interaction- induced effects. Let us note that the 
effect of beyond mean-field correlations is essential in 
this case. Mean-field calculations, either with realistic or 
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FIG. 12: (Color online) The = imaginary (left panels) and real (central panels) parts of E| as a function of the energy 
for different asymmetries at p = 0.16 fm~^ and T = 5 MeV. Upper (lower) panels display the results for neutrons (protons) 
obtained from the Avl8 interaction. The right panels display the partial derivatives of the real part with dots representing the 
location of the self-consistent quasi-particle energy. 



phenomenological forces, will predict results of the same 
order of magnitude of the FFG. Again, the latter pro- 
vides a measure of thermal effects. In the left panel, 
for instance, the correlated iso-depletion at the largest 
asymmetry considered should be significantly affected by 
thermal effects and therefore one should not take that 
value as a NN dynamical correlation effect. For the right 
panel, thermal effects are suppressed by degeneracy and 
the results can be taken as purely interaction-induced, 
zero-temperature values. 

Table I summarizes the values of n(0) for SNM and 
PNM at p = 0.16 fm-3 and T = 5 MeV for the different 
interactions considered in the left panel of Fig. 11. In all 
cases, n(0) for PNM is larger than for SNM, confirming 
the notion that, for the same density, neutron matter is 
less correlated than nuclear matter. The values of the de- 
pletion for both systems are within less than 5% of each 
other for the three phase-shift equivalent potentials (and 
also for Av6' and Av8'), which suggests that this quantity 
is rather well established at the theoretical level. Numer- 
ical uncertainties are well under control in this case and 
will only affect the third significant digit in the correlated 
calculations. Moreover, the difference of the occupation 
between neutron and nuclear matter is rather constant 
for all potentials. The conclusion is therefore warranted 
that the asymmetry dependence of the difference between 



the depletion for neutrons and protons is rather strongly 
constrained once an interaction is employed that yields 
an accurate description of the NN phase shifts. Av4' does 
not reproduce these and produces an unrealistically small 
depletion for SNM, thus effectively decreasing the differ- 
ence between nuclear and neutron matter. Due to the 
fact that both, neutron and nuclear matter, are degener- 
ate systems at this density and temperature, as indicated 
by the extremely small depletions for the FFG case, these 
results are not affected by thermal effects. As a matter of 
fact, the results obtained with Reid93 are in very good 
agreement with previous calculations at zero tempera- 
ture for the same potentials, obtained with a different 
procedure to implement the SCGF method [47]. 

The isospin asymmetry dependence of n^(0) can also 
be analyzed in terms of the time-ordered components of 
the self-energy. Figure 12 shows the fc = imaginary (left 
panels) and real (central panels) parts of S^, at p = 0.16 
fm~'^ and T = 5 MeV, as a function of the energy for 
different asymmetries. Upper (lower) panels display the 
results for neutrons (protons) obtained for the Avl8 in- 
teraction. The increase of isospin asymmetry has a repul- 
sive effect in the neutron ImE^, which is shifted by almost 
a constant amount as a increases. The minimum in this 
function becomes shallower and the threshold at lo ^ fXr 
is less pronounced. These repulsive effects are associated 
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with the decreasing importance of the tensor force for 
neutrons with increasing asymmetry. As a consequence 
of the diminishing threshold behavior, the inverted peak 
in ReS^ for neutrons is smeared out at large a's. The 
corresponding partial derivatives (right panels) become 
smaller in absolute value which, together with the quasi- 
particle shift, lead to a decrease of neutron depletion as 
asymmetry increases. Let us note here that the k = 
quasi-particle energies for neutrons become more repul- 
sive as asymmetry increases [58]. However, the distance 
with respect to the (also more repulsive) chemical po- 
tential increases with asymmetry, providing the shift to 
more negative values observed in the figure. 

In the case of protons, the situation for the down com- 
ponent is significantly different. The minimum of ImSj 
becomes larger as the system becomes more neutron rich 
and the threshold is more pronounced. Consequently, 
the inverted peak in the real part becomes larger and 
more attractive. When the partial derivative is taken 
(right panel), similar results are obtained for all asym- 
metries in the energy region of interest. This indicates 
that the effect of asymmetry on the proton self-energy 
is mainly energy independent, with the real part only 
shifting by a constant at energies oj < fir- This is in con- 
trast to the results for neutrons (upper left panels), where 
the partial derivative changes substantially as asymme- 
try increases, i.e. asymmetry induces energy-dependent 
effects for neutrons but not for protons. In spite of these 
qualitatively different behaviors, n„(0) and np{0) deviate 
from the symmetric value by almost the same amount as 
asymmetry increases (see Fig. 10). Similar results have 
been obtained for other NN interactions under the same 
conditions. In the view of these results, the "univer- 
sality" found for the iso-depletion for different interac- 
tions appears to be rather surprising. In spite of having 
very different ImSj^ components and predicting different 
quasi-particle shifts, the differences in the two momen- 
tum distributions are such that they give the same results 
independently of the NN force. 



VI. CONCLUSIONS 

The aim of the present paper is to provide a better 
understanding of the depletion of the nuclear Fermi sea 
as a function of density, temperature and isospin asym- 
metry, for different choices of realistic interactions. The 
SCGF method, when implemented at the level of lad- 
der diagrams, is capable of providing answers to these 
questions with emphasis on the role of short-range and 
tensor correlations that are induced by realistic NN in- 
teractions. While long-range correlations have an addi- 
tional influence on the depletion of the Fermi sea, it is ex- 
pected that they are more important near the Fermi en- 
ergy. For a proper comparison with finite nuclei, however, 
these long-range volume effects are presumably irrele- 
vant, since they must anyway be replaced by the surface 
dominated physics of low-lying nuclear states. The SCGF 



study reported here is performed at finite temperature, 
on the one hand to clarify the role of the temperature 
in determining depletion effects and, on the other hand, 
to avoid dealing with the technically challenging issue of 
pairing correlations. Since pairing is confined to energies 
very near the Fermi energy, at least at densities relevant 
for nuclei, our results will not be modified when such cor- 
relations will be considered in the future. We therefore 
report results mostly for T = 5 and 10 MeV, well above 
any pairing transition temperature for the systems un- 
der study [25]. We have studied both one-component 
systems (symmetric matter and pure neutron matter) as 
a function of density, as well as isospin-polarized matter 
as a function of asymmetry for normal and twice normal 
density. 

Our results indicate that temperature affects n(fc) in 
a qualitatively similar way as dynamical correlations do. 
The two types of correlations, however, induce quantita- 
tively different effects and can therefore be distinguished. 
To separate these components, we have relied on a com- 
parison with the FFG results, which are only affected 
by temperature effects. In the degenerate regime, i.e. 
at high density and relatively low temperatures, tem- 
perature effects are confined to a redistribution of parti- 
cles around the chemical potential within an energy scale 
that is proportional to the temperature considered. Dy- 
namical correlations induced by the short-range NN re- 
pulsion and tensor effects, on the contrary, generate an 
almost momentum-independent depletion of the Fermi 
sea, and, complementary, lead to an occupation of high- 
momentum states far from the Fermi momentum. 

The CDBonn and Avl8 interactions have been em- 
ployed to identify possible differences for the resulting 
nuclear momentum distributions. For symmetric nuclear 
matter and neutron matter, the depletion for fc = 
is strongly influenced by temperature at low densities. 
Once a density corresponding to p = 0.08 fm^'' is reached 
in PNM, short-range correlations lead to slowly increas- 
ing depletion as a function of density. This behavior is 
intuitive, since neutrons will be exposed more frequently 
to their mutual repulsion with increasing density, thus 
enhancing their occupation of high-momentum states, 
and, accordingly, increasing the depletion of the Fermi 
sea. For SNM, this behavior is not observed. The de- 
pletion saturates for the AvlB beyond normal density 
at a value of about 12%, while it continues to decrease 
for CDBorm. The larger depletion of the Avl8 as com- 
pared to CDBonn further emphasizes that the first is a 
slightly harder NN interaction. This is also illustrated 
by the different behavior of their self-energies. In or- 
der to understand this puzzling density dependence, we 
have utilized an approximate relation between the mo- 
mentum distribution and energy derivatives of the real 
part of the time-ordered self-energy taken at the corre- 
sponding quasi-particle energy (emphasizing the impor- 
tance of self-consistency). For momenta below kp, this 
result is valid when the energy dependence of the real 
part of the self-energy associated (through a dispersion 
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integral) with the imaginary part above the chemical po- 
tential is essentially linear. For momenta above kp, a 
similar relation holds when the real part of the self-energy 
is considered that is related to the imaginary part below 
the chemical potential. Near the Fermi momentum and 
at low densities, this approximation ceases to be valid. 
However, for most situations considered in this paper, it 
provides an excellent approximation to the momentum 
distribution calculated from the self-consistent spectral 
function. 

In order to utilize this relation, we have generated the 

time-ordered contributions to the nucleon self-energy at 
finite temperature. This allows us to gain access to the 
appropriate energy derivatives. Finite-temperature cal- 
culations are normally performed with retarded quanti- 
ties, but the "particle" and "hole" decomposition of the 
self-energy provides a more suitable connection with the 
corresponding results at zero temperature. The relevant 
self-energy components for determining the occupation 
inside the Fermi sea indeed shed light on the anomalous 
density dependence of the depletion in symmetric nu- 
clear matter. Mainly due to the tensor force, there is a 
an important threshold behavior in the imaginary part 
of the self-energy above the chemical potential in sym- 
metric matter, which is absent in neutron matter. For 
CDBonn, the threshold gives rise to an inverted peak in 
the corresponding real part that moves to higher energy 
with increasing density, yielding an almost density inde- 
pendent derivative. Since the tensor force also leads to 
more binding with increasing density, the location of the 
derivative moves farther away from the peak with density, 
thereby decreasing its absolute value and accordingly the 
depletion of the Fermi sea. For the Avl8 interaction, the 
energy derivative decreases with density on accoimt of its 
less dramatic threshold behavior for the imaginary part, 
yielding similar binding from an imaginary part that is 
significantly stronger than for CDBonn, but peaking at 
much higher energy. Since the binding effects are similar 
to those for CDBonn. the two effects compensate leading 
to the saturating behavior of the depletion observed in 
Fig. 2. For neutron matter, the ^Si-^Di coupled partial 
wave is absent and no strong threshold behavior is gen- 
erated in the imaginary part for both interactions. This 
leads to the more intuitive behavior that the depletion 
increases with the density. 

For asymmetric systems, we find a significant change 



of depletion with isospin asymmetry. We identify a tem- 
perature effect for the minority protons (but not for the 
majority neutrons) at an asymmetry beyond a = 0.4 at 
T = 5 MeV and normal density, due to the correspond- 
ing low density for protons. For twice normal density, 
the temperature effect does not play a role up to an ex- 
treme asymmetry of 0.8. The effect was implicit (and 
stronger) in earlier results obtained at T = 10 MeV [24] . 
The remaining difference between the neutron and pro- 
ton depletion is due to the decreased (increased) impor- 
tance of tensor correlations for neutrons (protons). By 
considering also different opcratorial clones of the AvlS 
interaction, like Av8', Av6', and Av4', we can indeed 
unambiguously demonstrate that this difference is asso- 
ciated with the tensor force. More importantly, this iso- 
depletion is independent of the chosen interaction, and 
therefore determined solely by phase shifts at least up to 
twice normal density. 

We close this discussion by noting that the employed 
interactions appear to slightly underestimate the deple- 
tion of the experimentally determined depletion of the 
deep proton mean-field orbits in ^°^Pb [16]. Interac- 
tions with an even stronger repulsion, like the Reid 
soft-core [73], yield a slightly larger depletion of about 
15% [27, 34, 74]. Nevertheless, while there remains a few 
% uncertainty as to the exact amount that is experimen- 
tally required, it is also clear that very different inter- 
actions, including modern ones, still lead to very simi- 
lar predictions for the depletion of the nuclear Fermi sea. 
Such depletions can be reliably calculated with the SCGF 
method. When lattice QCD calculations yield unambigu- 
ous information about the short-range NN repulsion, the 
remaining uncertainty about the depletion of the nuclear 
Fermi sea might eventually be eliminated. 
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